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ABSTRACT 

We study the effect of mergers on the morphology of galaxies by means of the 
simulated merger tree approach first proposed by Moster et al. This method combines 
N-body cosmological simulations and semi-analytic techniques to extract realistic ini¬ 
tial conditions for galaxy mergers. These are then evolved using high resolution hy- 
drodynamical simulations, which include dark matter, stars, cold gas in the disc and 
hot gas in the halo. We show that the satellite mass accretion is not as effective as 
previously thought, as there is substantial stellar stripping before the final merger. 
The fraction of stellar disc mass transferred to the bulge is quite low, even in the case 
of a major merger, mainly due to the dispersion of part of the stellar disc mass into 
the halo. We confirm the findings of Hopkins et ah, that a gas rich disc is able to 
survive major mergers more efficiently. The enhanced star formation associated with 
the merger is not localised to the bulge of galaxy, but a substantial fraction takes place 
in the disc too. The inclusion of the hot gas reservoir in the galaxy model contributes 
to reducing the efficiency of bulge formation. Overall, our findings suggest that merg¬ 
ers are not as efficient as previously thought in transforming discs into bulges. This 
possibly alleviates some of the tensions between observations of bulgeless galaxies and 
the hierarchical scenario for structure formation. 

Key words: galaxies: disc, evolution, interactions, structure - methods: numerical, 
N-body simulation 


1 INTRODUCTION 

Galaxy morphologies constitute one of the earliest attempts 
to classify galaxies, according to the relative prominence of 
their two main components, i.e. the (spheroidal) bulge and 
disc. This scheme was originally proposed 
. The nature of the link between the mor¬ 
phological properties of a galaxy and its cosmologically de¬ 
fined evolutionary track is a fundamental question in the 
field of galaxy formation and evolution. 

The standard model of galaxy formation assumes that 


* kannanr@mit.edu 
f maccio@mpia.de 


th e (exponential) 
bv lHubbld (Il926li 


the gas in dark matter potential wells cools to the center of 
the well and forms a disc, out of which stars form. Violent 
dissipative processes such as mergers and close encounters 
remove the angular momentum from the disc fuelling bulge 
formation. In a A cold dark matter universe (ACDM), dark 
matt er (DM) structures grow hierarchically (IWhite fc Reed 
[Hzi), with the smaller dark matter haloes forming first and 
later merging to form bigger ones: this makes interactions 
between galaxies a fundamental and inescapable process of 
galaxy evolution. Mergers are considered the origin of the 
so called “classical” bulges (i.e. whose properties are simi¬ 
lar to Elliptical gal axies, see e.g iDavies fc Illingworthlll983l : 
iDutton et al.ll2()l(^ . while “pseudo”-bulges (i.e. those char¬ 
acterized by disc-like profiles and/or kinematics, see e.g. 
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iKormendv fc Kennicuttll2004h are connected with “in-situ” 
processes like gravitational instabilities, which leads to the 
rearrangement of the disc material into a spheroidal-like 
structure. Additionally, a bulge-dominated galaxy may re¬ 
grow a new stellar disc, if there is a sufficient supply of 
cold gas, e.g., from the reservoir of hot gas present in quasi¬ 
hydrostatic equilibrium within th e gravitational poten tial 
of the dark matter halo (see e.g. iDe Lucia et al.l I 2 OIII ') or 
through cosmological accretion. The relative efficiency of 
these processes dictate the morphology of the galaxy. 

Many authors have tried to figure out the relative con¬ 
tribution of b ulge and disc components of galaxies in th e 
local universe llWeinzirl et al.l l2009l : iFisher fc DrorvI I 2 OIIII . 
iGadottil ll2009ll calculated the stellar mass content and dis¬ 
tribution for each galaxy component, and showed that, for 
galaxies more massive than 10^° Mq in the local universe, 32 
per cent of the total stellar mass is contained in ellipticals, 
and the corresponding values for discs, bulges and bars are 
36, 28 and 4% respectively. Classical bulges contain 25% of 
the total stellar mass, while pseudo-bulges contain 3% per 
cent. 

IKormendv et al] 1 I 2 OIOI) looked at galaxies in the local 
neighborhood (< 8 Mpc) and Hnd four galaxies consistent 
with being pure disc galaxies and 7 galaxies, including the 
Milky Way, having pseudo bulges. Given the dearth of clas¬ 
sical bulges in their sample, they thus estimate that around 
58 — 74% of the galaxies in their sample did not undergo vi¬ 
olent mergers in their past and thus they claim that the for¬ 
mation of these massive bulgeless galaxies represents a chal¬ 
lenge for current models of galaxy formation. A recent study 
of the m orphological mix in the S PSS volume has been dis¬ 
cussed in IWilman &: ErwinI ll2012h : they provide the fraction 
of galaxies showing a given morphological and activity classi- 
hcation as a function of stellar and parent halo masses. They 
Hnd that the fraction of elliptical galaxies is a strong func¬ 
tion of stellar mass; it is also a strong function of halo mass, 
but only for central galaxies. This is treated as evidence for 
a scenario where elliptical galaxies are always formed, prob¬ 
ably via mergers, as central galaxies within their halos, with 
satellite ellipticals being previously central galaxies accreted 
onto a larger halo. 

Different theoretical tools have been employed to ex¬ 
plain and understand this observational evidence. Bulge for¬ 
mation processes in Semi Analytic Models (SAMs) and their 
relative importa nce has been studied in detail in a number 
of recent papers llDe Lu cia et ai]l201Cll: i Fontanot et al.l[201ll : 


IWilman et al.l l2013l : IPorter et ai.ll2014l l. The general con¬ 

sensus in these works is that processes like mergers and disc 
instabilities are key to understanding the morphology of in¬ 
termediate mass galaxies (10'^° < M/Mq < 10^'^): therefore 
our limited understanding of these mechanisms is a limita¬ 
tion for the models’ ability of reproducing the morphological 
mix. In particular, IFontanot et al.l (1201 il l showed that the 
observed abundance of massive galaxies without a classical 
bulge is consistent with the predicted abundance of bulge¬ 
less galaxies only for a mo del where disc inst ability process 
is not considered (see also lPorter et al]l2014^ . 

If the angular momentum of a primordial halo is con¬ 
served during its collapse then it i s sufficient to produc e 
large discs (e.g., iFall fc Efstathioril Il980l : IMo et ahl Il998ll . 
Gooling in the highly dense inner regions and dynamical 
friction of orbital satellites dissipates the angular momen¬ 


tum of the gas (e.g., iD’Onghia et al.1 12 OO 6 I I , thus result¬ 
ing in compact discs in hydrodynamical N-body simula¬ 
tions, in which the rotation curve peaks at a few kpc, 
in contra diction to the rotation curves of observed gala x¬ 
ies (e.g., Steinmetz fc N avarro! Il999l: Maver et al. I 2 OO 8 I ) . 


ISpringel fc Hernauisd I 2005ll and iRobertson et alT ~( 2006fl 


showed that in idealized merger simulations with strong stel¬ 
lar feedback it is possible to get a disc dominated remnant, if 
the initial disc is gas rich. Furthermore, recent cosmological 


strong stellar 

feedback prescriotions llGovernato et al.1 

2007; 

Guedes et al. 

I 2 OIII: iBrook et al.l l2012l; Martig et al. 

2012; 

Stinson et al. 

2 OI 3 I: Kannan et al.1 20141: Marinacci et al. 


Christensen et al 


2OI4II . 


iHopkins et al.1 1 200E , hereafter H09) quantified the de¬ 
pendence of bulge formation on the gas disc fraction and 
presented a simple toy model to account for it. They showed 
that it is possible to obtain disc dominated remnants even 
for 1 : 1 gas rich mergers. This result has major implica¬ 
tions for the amou nt of bulge dominat ed galaxies found in 
the local universe dHopkins et al.ll2009fl . 

In this paper we re-examine the disc to bulge trans¬ 
formation dur ing mergers. We ado pt a hybrid method, first 
developed by iMoster et al.l (l2014l l which is based on high 
resolution hydrodynamical simulations of merger systems. 
The dark matter halo properties and their orbital param¬ 
eters are directly extracted from cosmological simulations, 
while the properties of the galaxies hosted by those halos 
are predicted using a Semi Analytic Model (SAM) of galaxy 
formation. 

In this way we gain the advantages of the merger simula¬ 
tions (high resolution and correct treatment of gas physics) 
and the SAM (cosmological background). Simultaneously, 
the computational cost is comparably low, so that a mean¬ 
ingful sample can be modelled in a short amount of time. 
Given the typical resolution of these hydrodynamical simu¬ 
lations, this method is well suited to resolve the small scales, 
relevant for the study of the evolution of the stellar compo¬ 
nents of galaxies and their scale parameters, such as the 
disc scale length and height, and for the evolution of galaxy 
morphology. This approach allows us to achieve the best res¬ 
olution possible within a reasonable amount of time, while 
being able to model a sample of galaxies in the correct cos¬ 
mological context. 

We describe the numerical techniques used to simulate 
mergers in (}2]a-iid also give a brief introduction of the SAM 
used in this study. We enumerate the results in 331 and the 
conclusions and discussions are given in fjd] 


2 MODELS 

In this section we briefly describe the methods that have 
been used in this paper. These are, the simulation code pino- 
CHIO that was employed to generate cosmological merger 
trees, the SAM morgana used to populate the merger 
trees with galaxies, the code to create initial conditions 
and the hydrodynamics code GADGET-2 used to perform the 
merger simulations. Throughout this paper we adopt cosmo- 
logical parameters c hosen to match results from WMAP-3 
(ISnergel et al.|[^007^ for a fiat ACDM cosmological model: 

= 0.26, Da = 0.74, h = 7/o/(100 km s"^ Mpc’^) = 0.72, 
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(Ts -— O.TT cLnd Ti — 0.95. Wg G.dopt Sr IlCrou 3 (l200lh IMF 
and compute all stellar masses accordingly. 


2.1 Merger Tree generation : PINOCCHIO 

To construct the Dark Matter halo merger trees we make use 
of the pinocchio code llMonaco et al.ll2002ll . pinocchio uses 
a scheme based on Lagrangian perturbation theory, which 
we b riefly describe in the following (see also iTaffoni et al.l 
I 2 OO 2 I for a more detailed discussion about the definition 
of DM haloes and merger trees). A Gaussian linear den¬ 
sity contrast field (for the desired cosmological background) 
is defined on a cubic grid, then smoothed repeatedly with 
Gaussian filters. For each particle on the grid, the six near¬ 
est particles are considered its Lagrangian neighbours. The 
collapse time of particles (i.e. the time they first belong 
to a high density, multi-stream region) is then computed 
for each point of the Lagrangian sp ace, following the de f- 
inition of orbit crossing proposed in [Monaco et al.l (l2002ll . 
Gollapsed particles accrete onto individual DM haloes or fil¬ 
aments according to their distance from the centre of mass 
of neighboring haloes (i.e. those containing at least one of its 
neighbours): if the distance is smaller than a given fraction 
of halo size, the collapsed particle then became an accreting 
particle of the halo. Filament particles can be accreted at 
later times if any of their neighbour became an accreting 
particle. Moreover, two DM haloes merge following a simi¬ 
lar prescription, i.e. if their distance in Lagrangian space is 
smaller than a fraction of the size of the more massive ob¬ 
ject. Overall PINOCCHIO allows a detailed reconstruction of 
the DM haloes, with known positions, velocities and angular 
momenta, and of their merger trees, in excellent a greement 
with the results of N-body simulations, (see e.g. Li et al.l 
1200111 . with a very fine time sampling that provides tracking 
of merging times without restriction to a fixed grid in time 
(as in N-body simulations). However, at variance to N-body 
trees pinocchio DM haloes are not allowed to decrease in 
mass, and the code does not track the evolution of DM sub¬ 
structures once they have been accreted by the main halo: 
however, these differences do not constitute a limitation in 
our case, since the mass evolution of DM substructures is 
explicitly tracked by the hydrodynamical simulation. 


2.2 Semi-analytic model: MORGANA 

The PINOCCHIO merger trees have then been used as 
input for the Sem i Analytic Model (SAM) MORGANA 
llMonaco et al.l[2007l L In SAMs, the evolution of the bary- 
onic component is followed by means of approximate, yet 
physically grounded, analytic prescriptions for modelling the 
relevant processes (such as gas cooling, star formation and 
feedback) and their interplay, as a function of the physical 
properties of model galaxies (like their cold gas and stel¬ 
lar content). These analytic prescriptions involve a number 
of parameters, usually calibrated by comparing model pre¬ 
dictions with a well defined set of low-redshift observations. 
Despite (and thanks to) this simplified approach, SAMs have 
turned into a flexible and powerful tool to explore a broad 
range of specific physical assumptions, over scales that could 
not be directly modelled simultaneously (ranging from the 
accretion onto a super-massive black hole on sub-pc scales 


to the Mpc scales involved in cosmological structure forma¬ 
tion) . In the following we will briefly describe the treatment 
of the most rel evant processes leadin g to bulge formation in 
MORGANA fsee lPe Lucia et al.ll201l] . for a discussion of the 
different channels for bulge formation in different SAMs). 

MORGANA distinguishes between minor and major 
galaxy mergers; the threshold of the two events being defined 
by a mass ratio between secondary and primary galaxy of 
0.3. The orbital decay of dark matter subhaloes and galaxy 
mergers are modell ed using the fitting formulae defined by 
ITaffoni et al.l (l2003l i. In case of a minor merger, the stellar 
mass and the cold gas of the secondary galaxy are completely 
given to the bulge component of the remnant galaxy, while 
the disc is considered unaffected. On the other hand, in case 
of major mergers, the whole stellar and gaseous disc of the 
merging galaxies are destroyed and relaxed into a spheroidal 
remnant. In both cases, any cold gas eventually associated 
with the bulge can be efficiently converted into stars, on very 
short time-scales (effectively triggering a starburst). At later 
times, the remnant spheroidal galaxy can grow a new disc, if 
cooling processes are effective in the parent halo. In partic¬ 
ular for this work, we make use of the standa rd realization 
of MORGANA defined in iDe Lucia et aP ll201lli . 


2.3 Galaxy models for N-body simulations 

We use the method described in ISoringel et al.l (l2005[l to 
initialize our galaxies. Each object is composed of five dif¬ 
ferent components: (i) a cold gaseous disc with mass Meg, 
(ii) a stellar disc (Maisc ), (in) a stellar bulge (Mb ), (iv) a 
gaseous halo (Mhg ) and (v) a dark matter halo (Mdm ). 

The gaseous and stellar disc have an exponential surface 
density profiles and their scale lengths (rg and ra respec¬ 
tively) are related via rg = yra, with x = l-S- We adopt a 
sech^ profile with a scale height zo for the vertical structure 
of the stellar disc and the vertical velocity dispersion is set 
equal to the radial velocity dispersion. The vertical struc¬ 
ture of the gaseous disc is computed by requiring balance 
between the galactic potential and the pressure given by the 
Equation of State (EOS); the EOS also fixes the tempera¬ 
ture of the gas, rather then t he velocity dispe rsion. Finally 
we adopt an Hernquist profile iHernauisli (Il990fl for both the 
stellar bulge and the dark matter halo. The stellar bulge has 
scale length rb, while the dark matter halo is defined by a 
scale length Ts, a concentration parameter c = Tvir/ra and a 
halo spin A. 

To model the hot gaseou s hale we f ollowe d 
the same parameterization of iMoster et al.l ll201l1) . 
Namely we use the obse r vatio n al motivated /3-profile 
l [Cayaliere fc_Fusco-Femiand Il976l . Ijones fc Formaiil Il984l . 
lEke et ah 1998h : 


phg(r) = po 



( 1 ) 


We use 0 = 2/3 lljones fc FormanI ll98J '). Tc = 0.22rs 
l|Makino et al.l[T998l i and fix po such that the hot gas mass 
within the virial radius is Mhg. The temperature profile is 
fixed by imposing hydrostatic equilibrium. In addition, we 
impose a slow rotation around around the spin axis of the 
discs for the hot halo. The specific angular momentum is 
defined as jhg and is set to be a multiple of the specific 
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angular momentum of the dark matter halo jdm such that 
Jhg = rrjdm- 

High _ reso lution cosmological simulations 

llGovernato et all |2010| . e.g.) have indicated that this 
“spin factor” a is generally larger than unity, as feedback 
processes preferentially remove low angular momentum 
material f rom the halo. For the exact value we have used the 
results of iMoster et alJ (l 201 ll '). who obtained a = 4 using 
isolated simulations of a MW-like galaxy and requiring that 
the evolution of the average stellar mass and scale-length 
found observationally be reproduced. 

2.4 Simulations 

The hydrodynamical simulations have been p erformed 
with the parallel TreeSPH-code GADGET-2 (ISoringell 
l2005h. The code uses Smoothed Partic le Hydrodynamic s 
I Lucvlll97^ ; iGingold &: MonaghanlllQT^ ; lMonaghSll992l ~l 
to evolve the gas using a n entropy conserving scheme 
dSpringel fc Hernauisdl2002l l. The code includes Radiative 
cooling f or a primor d ial m ixture of hydrogen and helium 
following iKatz et al.l (Il996l l and a spatially uniform time - 
independent local UV background jHaardt fc Madaulll996ll . 
We model Star formation and the associated heating by 
supernovae (SN) follo wing the sub-resolution mul tiphase 
ISM model described in ISpringel fc Hernauisd (120031 '). Cold 
clouds form stars in dense (p > Pth) regions on a time scale 
chosen to match observations (iKennicu ^ Il998l ') . The thresh¬ 
old density pth is determined self-consistently by demanding 
that the equation of state (EOS) is continuous at the onset 
of star formation. 

We also include SN-d riven galactic winds following 
ISpringel fc Hernauisd (l2003f ). The wind mass-loss rate is as¬ 
sumed to be proportional to the star formation rate (SFR) 
Mu, = CM*, where the mass-loading-factor C quantifies the 
wind efhciency. We assume a constant wind speed model 
with Vw = 500kms~ (energy-driven wind). We refer to 
Table [1] for the values of further parameters (assuming a 
Kroupa IMF) describing the multiphase feedback model. We 
do not include feedback from accreting black holes (AGN 
feedback) in our simulations. 

2.5 Simulations of Semi-Analytic Merger Trees 

We use the method devised bv iMoster et al.l |(20l3) to com¬ 
bine (pinocchio)-I-SAMS results with high resolution hy¬ 
drodynamical simulations. As a first step we select a merger 
tree from the pinocchio simulation, we then use the MOR¬ 
GANA semi-analytical model to predict the properties of the 
galaxies hosted in each branch of the tree. The prediction 
for the properties of the galaxies directly taken from the 
MORGANA SAM include the masses of the gas disc, stel¬ 
lar disc, stellar bulge, hot gaseous halo and the dark matter 
halo and the scale lengths of the stellar disc and the stellar 
bulge. These predictions are used to create a particle based 
realization of the galaxy (made of dark matter, stars and hot 
and cold gas, see section [?!3ll at the designed started time of 
the simulation (here Zi = 1). The number of particles in each 
component is decided by fixing A* = 500, 000, which is the 
total number of stellar particles in the final merger remnant. 
This then sets the mass of the star particles. We also impose 


that the mass of DM particle be 15 times more massive than 
the stars and the gas particles to be 2 times more massive. 
The typical mass of DM, gas and stars in our simulations are 
1.0 X 10® M 0 , 1.4 X 10® Mq and 6.6 x 10"* Mq respectively. 
The softening length (ci) of each component ‘i’ of the galaxy 
is given by = 32 kpc x y/mpart.i/10i° M©, where mpart,i 
is the mass the particles of the ‘i’th component. Typically for 
our simulations this equates to a softening of about about 
100 pc for the gas particles, 300 pc for the DM particles and 
80 pc for the stars. 

This galaxy is then evolved with the hydrodynamical 
code until the time of its first merger, as predicted by the 
merger tree. At this point we stop the hydro simulations and 
create a particle realization of the satellite system, which is 
also based on the semi-analytic prediction. We then add the 
satellite in the simulation at the virial radius of the main 
halo and we restart the hydrodynamical run. The orbital 
parameters (position and velocity at the time of accretion) 
are directly taken from the A-body simulation, in this way 
“naturally” creating a cosmologically motivated merger. The 
system galaxy-l-satellite is evolved with the hydrodynamical 
code, until the next satellite galaxy enters the main halo. 
This process is repeated for all merging satellites until z = 0. 
Table [1] lists all parameters used to construct and run our 
galaxies (see also lMoster et all (l2014fl for further details and 
the exact meaning of all parameters). 


3 MORPHOLOGICAL EVOLUTION OF 

GALAXIES 

As mentioned earlier the most important processes which 
dictate the morphology of galaxies are mergers and close en¬ 
counters. These processes are known to trigger mass transfer 
into the central bulge through different channels; the main 
ones being : (i) accretion of satellite material onto the bulge 
of the central galaxy (central bulge) during a merger, (ii) the 
transfer of stars from the disc of the central galaxy (central 
disc) to central bulge and (iii) funnelling of gas towards the 
centre and subsequent star formation. 

As described in Sec. SAMs use simple prescriptions 
for mass transfer through these channels. Are these pre¬ 
scriptions correct? For example the simple assumption that 
during a major merger all the material of t he central disc 
goes i nto the bulge has been contradicted bv iHopkins et all 
ll2009ll . who showe d that discs can sur vive major mergers 
if they are gas rich. IChang et al.l ll2013ll showed that a disc 
dominated satellite is easily tidally stripped as it orbits a 
halo before it finally mergers, which reduces the amount of 
material given by the satellite to the central bulge. 

Here we test different scenarios using high resolution hy¬ 
drodynamical simulations of galaxy mergers, starting from 
the initial conditions given by the semi-analytic model as 
described in the previous section. We simulate a total of 20 
merger events covering a range of different merger histories 
and galaxy properties. The merger trees and their parame¬ 
ters are given in Table [2] Finally for each tree we also per¬ 
form an ’isolation’ run, i.e. a run where the central galaxy 
is evolved without any mergers from z = 1 to the present 
time. The isolated runs act as control, which helps to easily 
dis-entangle the effect of mergers. 
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Table 1. Summary of the parameters used for the simulations of merger trees and their fiducial value. 


Parameter 

Description 

Fiducial value 


Redshift at the start of the simulation 

1.0 

Mmin 

Minimum dark matter mass ratio 

0.03 

<5 

Ratio of scaleheight and scalelength of the stellar disc 

0.15 

X 

Ratio of scalelengths between gaseous and stellar disc 

1.5 

? 

Ratio of gaseous halo core radius and dark matter halo scale radius 

0.22 

/3hg 

Slope parameter of gaseous halo 

0.67 

a 

Ratio of specific angular momentum between gaseous and dark halo 

4.0 

N. 

Expected final number of stellar particles in the central galaxy 

200 000 

K. 

Ratio of dark matter and stellar particle mass 

15.0 

-^reSjSat 

Ratio of satellite and central galaxy particle mass 

1.0 

N ■ 

^ ^min 

Minimum number of particles in one component 

100 

ei 

Softening length in kpc for particle of mass m = 10^*^ Mq 

32.0 

i* 

Gas consumption time-scale in Gyr for star formation model 

3.5t 

^0 

Cloud evaporation parameter for star formation model 

1250.0+ 

Psv 

Mass fraction of massive stars for star formation model 

0.16+ 

7s N 

Effective supernova temperature in K for feedback model 

1.25 X 10®+ 

C 

Mass loading factor for wind model 

1.0 

rJ-wind 

Initial wind velocity in kms“^ for wind model 

500.0 


t The star formation parameters assume a Kroupa IMF. 


3.1 Bulge/Disc Decomposition 

There are many different ways to decompose the mass of 
a galaxy into its basic morphological entities i.e., disc and 
bulge, for example by fitting surface brightness profiles wit h 
a bulge and disc component (e.g., iRobertson et al.ll2006ll . 
kinematic decompositions based on one or two dimensional 
kinematic maps and three dimensional component fits. All 
these methods rely heavily on the viewing angle. In this 
paper we make use of the six dimensional phase space of 
(x,y,z,Vx,Vy,Vz), to track the bulge and the disc compo- 
nent s of the primary galaxy, t hroughout the simulatio n (see 
also IScannapieco et Si l2010l : iMarinacci et al.l l2014lJl . The 
galaxy in question is viewed edge on i.e., the angular mo¬ 
mentum vector of the disc is placed parallel to the 2 : axis. 
Now for a purely rotating disc the circular velocity of a given 
particle, at a distance r is given by 


GM{< r) 


( 2 ) 


and the specific angular momentum of the particle will be 


lc = r 


GM{< r) 


(3) 


If a particle is purely rotationally supported then the 
ratio of its specific angular momentum in the z direction 
(G) and G if rot) will be equal to 1. 


frot — 


(4) 


If we consider all the stellar component of a sample 
galaxy from our simulations and plot the distribution of the 
mass as a function of the rotational support i.e., G/G (Fig- 
1(a) I, we expect it to be bi-modal with the rotationally sup¬ 
ported disc particles distributed around frot ~ 1 and the 
velocity dispersion (a) supported bulge particles grouping 
around frot ~ 0. Decomposing the distribution into bulge 


and disc components comprises of figuring out the local max¬ 


ima close frot = 0. In the simple example shown in Fig. 1(a) 


frot is symmetrical around 0, meaning that the bulge is non¬ 
rotating. This might not be true for many galaxies where the 
bulge might have some low amount of rotation induced dur¬ 
ing a merger (for e.g. see right panel of Fig. 2(b) I. In order 
to make a self consistent decomposition, instead of assum¬ 
ing the symmetric nature of frot around 0, we calculate the 
point at which the slope or first derivative of the stellar dis¬ 
tribution function becomes zero, particularly points where 
the slope changes from positive to negative (i.e. local max¬ 
ima). This takes care of cases where the bulge has a rotation, 
but this only works if the bulge and disc distributions are 
quite separated from each other. In some galaxies the two 
distributions merge and the function no longer has a local 
maxima. This problem is overcome by calculating the in¬ 
flection points of this function, especially the points where 
the second derivate changes sign from negative to positive 
(i.e. the curve changes from convex to concave). This neatly 
separates the bulge and disc components even in galaxies 
where the stellar and bulge distributions overlap quite a bit. 
It should be mentioned that we look for these local max¬ 
ima and inflection points only between —0.3 ^ frot G 0.3, 
beyond which we consider that the galaxy has no bulge 
component. The stars in the bulge are then assumed to be 
symmetrically distributed around the local maximum or in¬ 
flection point, while the disc stars make up the rest of the 
distribution (as shown in Fig. |l(a)1 l. This technique, by con¬ 
struction has the ability to only statistically determine the 
bulge and disc components of a simulated galaxy and cannot 
uniquely classify individual particles as belonging to either 
of the components. For example in Fig. 1(a) all particles 
with frot ~ 0.7 have equal probability of being a bulge and 
disc particle. 


This decomposition on a sample galaxy at redshift z=l, 
from our catalog (Tree 350) gives a B/T ~ 0.7, which is 
confirmed by the initial conditions taken from the SAM. 
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Table 2. Table listing the properties of simulated merger trees 


Tree/Sat ID 

^enter^ 



log{Mh)^ 

log(NU)^ 

log{Mcg]l 

r£ 

Tree 18989 

1.0 

- 

- 

11.70 

9.92 

10.06 

- 

Sat 1 

0.98 

0.37 

0.83 

11.27 

9.92 

9.65 

0.20 

Tree 28678 

1.0 

- 

- 

11.67 

10.44 

9.62 

- 

Sat 1 

0.46 

0.56 

0.16 

11.47 

9.66 

8.74 

0.44 

Tree 80891 

1.0 

- 

- 

11.86 

10.61 

9.65 

- 

Sat 1 

0.81 

0.05 

0.037 

10.67 

8.80 

8.92 

0.81 

Tree 65521 

1.0 

- 

- 

11.81 

10.48 

10.21 

- 

Sat 1 

0.77 

0.10 

0.05 

10.84 

9.05 

8.93 

0.43 

Tree 154448 

1.0 

- 

- 

11.631 

10.60 

9.83 

- 

Sat 1 

0.67 

0.97 

0.74 

11.651 

10.74 

8.79 

0.42 

Tree 215240 

1.0 

- 

- 

11.54 

9.84 

9.88 

- 

Sat 1 

0.37 

0.76 

0.47 

11.52 

9.86 

8.21 

0.13 

Tree 455141 

1.0 

- 

- 

11.39 

10.14 

9.54 

- 

Sat 1 

0.56 

0.88 

0.38 

11.45 

9.48 

9.54 

0.56 

Tree 114590 

1.0 

- 

- 

11.15 

9.16 

8.96 

- 

Sat 1 

0.92 

0.86 

0.84 

11.10 

9.12 

9.32 

0.71 

Tree 28837 

1.0 

- 

- 

11.86 

10.42 

9.94 

- 

Sat 1 

0.72 

0.04 

0.01 

10.52 

8.48 

7.33 

0.31 

Tree 50967 

1.0 

- 

- 

11.86 

10.34 

10.17 

- 

Sat 1 

0.45 

0.062 

0.05 

10.69 

8.46 

8.52 

0.43 

Tree 52201 

1.0 

- 

- 

11.64 

10.01 

10.03 

- 

Sat 1 

0.34 

0.66 

0.51 

11.53 

9.76 

8.825 

0.75 

Tree 53334 

1.0 

- 

- 

11.78 

10.23 

9.58 

- 

Sat 1 

0.72 

0.45 

1.0 

11.47 

10.22 

9.38 

0.82 

Tree 58811 

1.0 

- 

- 

11.78 

10.36 

9.65 

- 

Sat 1 

0.17 

0.38 

0.3 

11.40 

9.55 

9.49 

0.15 

Tree 60367 

1.0 

- 

- 

11.83 

10.47 

9.98 

- 

Sat 1 

0.39 

0.03 

0.007 

10.41 

8.40 

7.42 

0.77 

Tree 61557 

1.0 

- 

- 

11.78 

10.40 

10.00 

- 

Sat 1 

0.28 

0.34 

0.13 

11.37 

9.53 

8.86 

0.14 

Tree 102663 

1.0 

- 

- 

11.92 

10.76 

10.10 

- 

Sat 1 

0.79 

0.27 

0.09 

11.37 

9.06 

9.62 

0.51 

Tree 350 

1.0 

- 

- 

11.90 

10.77 

9.84 

- 

Sat 1 

0.96 

0.13 

0.041 

11.02 

9.25 

8.67 

0.35 

Sat 2 

0.95 

0.04 

0.016 

10.53 

8.75 

8.52 

0.21 

Sat 3 

0.93 

0.04 

0.009 

10.62 

8.67 

7.73 

0.36 

Tree 2536 

1.0 

- 

- 

11.63 

10.36 

10.15 

- 

Sat 1 

0.79 

0.12 

0.025 

10.75 

8.76 

8.28 

0.69 

Sat 2 

0.75 

0.29 

0.312 

11.20 

9.98 

6.20 

0.66 

Sat 3 

0.46 

0.19 

0.135 

11.18 

9.46 

8.94 

0.31 

Tree 187460 

1.0 

- 

- 

10.94 

9.14 

9.21 

- 

Sat 1 

0.70 

0.80 

0.85 

11.09 

9.17 

8.82 

1.05 

Sat 2 

0.61 

0.66 

0.79 

11.32 

9.26 

8.10 

0.63 

Tree 159419 

1.0 

- 

- 

11.05 

9.08 

8.26 

- 

Sat 1 

0.89 

0.52 

0.41 

10.79 

8.75 

7.80 

0.18 

Sat 2 

0.88 

0.77 

0.66 

11.15 

9.34 

8.15 

0.49 


^ Redshift of entry of the halo 
^ Dark Matter merger ratio 
^ Baryon merger ratio 
^ Dark Matter mass 
^ Stellar mass 
f Cold gas mass 
9 Circularity parameter 


The stellar density maps (figure 1(b) I also confirm that this 
is a bulge dominated galaxy. 


Fig. 2(b) shows the effect of a minor and major merger 
events and how the angular momentum distribution of mat¬ 
ter in the central galaxy changes in our simulations. A mi¬ 
nor merger does not particularly affect the morphology of 


the galaxy (Fig 2(a) I which remains disc dominated (Tree 


65521). On the other hand a multiple major merger com¬ 
pletely changes the galaxy morphology as shown in of Fig. 
2(b) A disc dominated galaxy turns into bulge dominated 


one with a very small hint of a remnant disc (Tree 187460). 
This method also allows us the track the mass flow into 
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Figure 1. An sample decomposition of a bulge dominated galaxy in our simulations. 

a) The mass distribution of the stellar component of the galaxy as a function of the rotational support. There is clearly a bi-modal 
distribution, with clear difference between the bulge (green) and disc (red) components. This particular decomposition yields a B/T ^ 0.7. 

b) The column density of the stellar component shown in projection, perpendicular and parallel to the total angular momentum axis of 
the galaxy. 


the central bulge during the whole simulation, as shown in 
Fig. El (Tree 65521; /r = 0.1) and|l(Tree 18989; /r = 0.76). 
Every particle in a GADGET-2 simulation has a unique ID, 
and for DM and stars, this ID does not change as a function 
of simulation time. When we set up the initial conditions 
we decompose the galaxy into bulge/disc components and 
make a note of the IDs of the particle in each component. 
As the simulation evolves in time, the particles are tracked 
through their unique IDs as they move from one component 
to other. This allows us to track the various mass transfer 
channels very efficiently. Since we only look at particles that 
are present at the start of our simulation we do not track 
how the morphology of the newly formed stars changes. 

We track the stars of the satellite in the same way and 
we consider that the particle has been unbounded from the 
satellite if the particle is at a distance > 3rs, or lies beyond 
5 kpc above or below the disc. Some extreme galaxies might 
have a massive bulge and might have bulge stars above the 
5 kpc limit that we have used here. In order to see if this 
classihcation is robust we looked at satellite stripping with 
the stripping radius of 7.5&10.0 kpc above and below the 
disc of the galaxy. In all the simulations of our sample we 
get less than 10% changes in all related merging channels 
when we make this change. This is because even if there is 
some amount of bulge stars present beyond this height, the 
density of stars is very low. This makes the error on the 
channels also quite low. Therefore, we continue to use this 
definition throughout the paper. 

Figures [3] & |4] show the evolution of various properties 
as a function of simulation time. The top left panel in these 


hgures show the evolution of the bulge/total ratio of the 
central galaxy in both merger and isolation. During a mi¬ 
nor merger the difference between the B/T ratio is minimal 
whereas when a galaxy undergoes a major merger the B/T 
ratio increases considerably. The middle panel on the left 
shows the distance of the satellite to the central galaxy. Due 
to the dependence of the merging time on the mass ratio of 
the merger, a low mass satellite spends more time orbiting 
the primary galaxy. The bottom panel on the left shows the 
mass loss of the satellite. The low mass satellite undergoes 
numerous close encounters, losing most of its mass to the 
central halo in the process, whereas the high mass satellite 
has fewer orbits and deposits almost all of its mass into the 
central bulge. The top panel on the right shows the amount 
of satellite stellar mass given to the central bulge and it is 
seen that a minor satellite deposits a very small fraction of 
its stellar mass into the central bulge owing to the fact that 
most of its mass is stripped away and lost to the halo during 
its descent while a massive satellite deposits almost all of its 
mass into the central bulge as it spends less time orbiting 
and it is also able to hold on to more of its mass before it 
merges. The middle panel on the right shows the evolution 
of the amount of central stellar disc mass given to the cen¬ 
tral bulge during a merger and in isolation. There is almost 
no mass transfer from the central disc to the central bulge 
in case of a minor merger, but during a major merger about 
30% of the central stellar disc ends up in the central bulge. 
This number is very small compared to most SAM prescrip¬ 
tions which assume that almost all the stellar mass of the 
central disc ends up in the bulge during a major merger. 
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Figure 2. The angular momentum distribution of the central galaxy before(blue curve) and after(green curve) a merger event. Left 
Panel: Change in stellar mass distribution of the central galaxy after undergoing a 1:10 merger. The merger heats up the thin disc, 
leaving the bulge almost completely unaltered. Right Panel: Change in stellar mass distribution of the central galaxy after undergoing 
two major mergers (1:1.25 & 1:1.5). Due to the high merger ratios, the galaxy builds up a large bulge, with only a very small hint of a 
remnant disc. 


As we will discuss later a large fraction of the central disc 
gets dispersed into the stellar halo, and a small remnant disc 
survives. Finally, the bottom panel on the right shows the 
composition of the central bulge. As expected the low mass 
satellite contributes little, while a high mass satellite con¬ 
tributes a large amount of material to the bulge and, in the 
end, it dominates the bulge composition in this particular 
example. 

In the following sections we look closer to the various 
channels for mass transfer into the bulge and we will also 
try to quantify them with simple, physically motivated, em¬ 
pirical relations. 

3.2 Importance of merger parameters 

The major parameters which control the outcome of a 
merger event is the ratio of masses of a merging galaxies 
and the their ration of orbital angular momentum. The ra¬ 
tio of the total halo mass of the galaxies (jj.) will control the 
dynamics of infall, while the components which dominate 
the potential at the center of the halos (stars -I- cold gas) 
will provide most of the torque during the final stages of the 
merger. Hence, we can expect the baryonic (stars -I- cold 
gas) merger ratio (fib) to be important as well. In addition 
we also the dependence of various transfer channels on fic, 
the ratio of the amount of total material within two scale 
radii (calculated here assuming that the DM halo follows a 
NFW profile) of the halos. This parameter is chosen so as 
to compare our results to H09. _ 

In addition, previous studies dChang et al.l[20l3l have 
stressed the importance of the satellite orbital parameters 
in determining the final outcome of the merger. In our ap¬ 
proach we extract our merger trees and the satellite orbital 
parameters directly from cosmological simulations. This en¬ 
sures that the orbital properties of satellites fall within the 
gamut of the normally occurring merging paths in the Uni¬ 
verse. 

The drawback of this method is that, there is no con¬ 


trol over the input parameters. This is makes it difficult 
to make a controlled experiment on the impact of vari¬ 
ous merger parameters (see fo r example the CALMER data 
base: [Chilingarian et al.ll2010l ). However it is still important 
to check the validity of previous works in a full cosmological 
setting as we do. However, we do try to quantify the effects 
of the mergers on the galaxy morphology as a function of 
various parameters. Our analysis is not meant to fully repro¬ 
duce the dependence of mass transfer on the merger prop¬ 
erties (our sample is too small for such a task), but to grasp 
the main trends and dependencies and provide guidance to 
test the implications of our findings in broader context. 

We look at the most important factors that could af¬ 
fect the mergers. A high mass satellite experiences a high 
dynamical frictional force and will fall quickly to the center, 
thereby depositing most of its mass into the bulge of the 
central galaxy. The opposite is true for low mass mergers. 
The satellite orbits around the central galaxy many times, 
losing its mass through repeated close encounters before it 
merges. The circularity of the orbit can be expressed as a 
function of the ratio between the orbital angular momen¬ 
tum of the satellite (j) and jc{E) the angular momentum of 
a circular orbit with the same energy as the satellite, i.e.. 


Of course the mass transfer in various channels will be di¬ 
rectly proportional to the merger ratios fi, Hb and/or 
and inversely proportional to rj. We quantify the dependence 
on orbit by exp(1.9rj), where the exponential function and 
the proportionality factor of 1.9 comes from the relation be- 
tween n and the total merg er time for satellites derived by 
iBovlan-Kolchin et al.l (120081 ) from cosmological simulations. 
We test the dependence of mass transfer through various 
channels on the parameters. We use a very simple functional 
form 
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Figure 3. Time evolution of merger parameters for Tree 65521 (/i = 0.10). Panel 1 - Evolution of B/T in both the isolated (dotted) and 
merger (solid) runs. Panel 2 - Distance of the satellite form the centre of the primary galaxy, Panel 3 - Mass loss of the stellar component 
of the satellite as it orbits in the primary halo of the central galaxy. Panel 4 - Stellar mass of the satellite given to bulge of the primary 
galaxy, Panel 5 - Amount of initial disc mass given to the bulge, (during a merger: solid curve and in isolation: dotted curve) and Panel 6 
- Composition of the bulge. The brown line denotes the fraction of stars which existed in the central bulge at the start of the simulation, 
the blue line is the contribution from the satellite and the black (solid and dotted) lines are the contribution from the central disc during 
a merger and in isolation. 


f{x) = ax 

where x e {/i, Ub, Uc, 


Uh 


f^C 


ea;p( 1 . 9 r 7 )’ 62:^(1.977)’ 62:^(1.977) 


} 

( 6 ) 

which ensures that the boundary value condition of 


/(x) = 0 at X = 0 is satisfied. For each of the mass transfer 
channels we test its dependence on all the four parameters 
and six functional forms. Table [3] shows the best fit value 
of a for all the parameters and it also gives the correspond¬ 
ing . We see that the results are not very different for the 
different merger parameters; moreover, a more sophisticated 
analysis designed to differentiate between the parameters is 
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Figure 4. Time evolution of merger parameters for Tree 18989 (p = 0.37). The panels denote the same quantities as in Figure|3] 


not feasible given our limited sample. In fact, the analysis 
doesn’t provide compelling evidence in favour of adding pa¬ 
rameters and/or physical dependencies, with respect to the 
simple assumption of 1-parameter linear dependence from a 
given physical property. Therefore, we just consider depen¬ 
dencies from the parameter providing the overall smaller x^ 
value, which in our case happens to be fib , and plot the 
functional dependence of the transfer channels on fib from 
now on. 


3.3 Where does the satellite mass end up? 

One of the major channels to build up the bulge of a galaxy 
is through the transfer of material from the satellite directly 
into the bulge. The morgana model employs a simple set of 
formulae to evolve the bulge mass during mergers. For ex¬ 
ample it is assumed that the whole satellite mass is added to 
the bulge, irrespective of the mass ratio between the satellite 
and the host and/or the orbit of the satellite. 

We define fab as the fraction of the initial (before infall) 
stellar mass of the satellite given to the central bulge. Fig¬ 
ure [5] show the results for fsb from our suite of simulations 
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Figure 5. The fraction of satellite mass given to the central bulge 
as a function of their baryonic merger ratios. The green curve is 
out fit to the simulated data (blue stars). 

(blue stars) as a function of ^t,. The simulations data show a 
significant dependence of fst on the baryonic merger ratio. 

The best fit parameters that describe fsb in our simu¬ 
lations are given by: 

fsb = 1.0/ii, (7) 

represented by the green line in figure O which is a fair rep¬ 
resentation of the simulations results (also matches results 
from H09), while being still simple to implement in semi- 
analytic models. 

3.4 Where does the central disc mass end up? 

Another major channel to build up the bulge is the transfer 
of mass (stars-hgas) from the central disc into the central 
bulge during and after a merger. In the MORGANA code the 
following formula is used to determine the mass (Mdb) trans¬ 
ferred through this channel: 

Mdb f 0 p ^ 0.3 

^ I 1 M > 0.3 ^ ^ 

where Mdiac is the central disc mass (stars and gas). On the 

other hand the gaseous and stellar components of the disc 
behave differently as pointed out by H09. They have shown 
that if the central disc is gas rich then only a small amount 
of its initial gas disc is funnelled into the bulge. This mass 
transfer is given by the following relation 

fcgb = (1 “ fgas)fJ'C (9) 

where fcgb is the fraction of the central disc gas mass trans¬ 
ferred to the central bulge, and fgas is the gas fraction of the 
central disc (ratio between gaseous and stellar disc mass). 

is the central (within two scale radii) merger ratio as de¬ 
fined by H09 and used in man y S AM calculations (s ee for 
example ISomerville et al]l2012l and IPorter et al.]l201^ . 

In this section we show the results of single merger 
events in our simulations inventory. The effect of multiple 
mergers on the central disc is very hard to quantify. The 
mass transfer from the disc to bulge, seems to be dependent 



Figure 6. The fraction of central gas disc mass given to the 
central bulge. The blue points are simulation results and the green 
line is the empirical fit to the points. 

not only on the properties of the merging galaxies but also 
on the time lag between the mergers because it is easier to 
perturb an already unstable disc. Therefore, quantifying the 
effect of multiple mergers is quite involved and we defer this 
to a future paper. In this work we only concentrate on the 
results from the binary merger simulations, the results of 
which are shown in Fig.[6l We find an empirical relation for 
the amount of gas mass transfer into the bulge given by 

fcgb = (1 “ fgas)Ub (10) 

shown as solid (green) line in figure [G] Therefore, we con¬ 
firm previous results on the importance of the gas fraction 
in determining the gas transfer from central disc to central 
bulge. 

On the other hand we find quite a substantial difference 
with previous studies for what concerns the fate of the cen¬ 
tral stellar disc after a merger. Many SAMs assume that all 
the stellar mass of the central disc ends up in the central 
bulge after a major merger, while H09 found that this mass 
transfer is given by 

fdb = Uc (11) 

However in our simulations even for a 1 : 1 merger the 
amount of central stellar disc mass that ends up in the cen¬ 
tral bulge is only about 37%, as shown in Fig. [7] The em¬ 
pirical fit to simulation data points is gives us a relation 

fdb = 0.37Mb (12) 

where fdb is the fraction of central stellar disc mass that 
ends up into the bulge. 

This result is very interesting and diverges a lot with 
the simplest picture of major mergers destroying discs, in¬ 
cluded in most (but not all) SAMs. Now the question arises: 
where does the rest of the disc end up? Major mergers are 
very violent events, where a lot of energy (mainly orbital) 
is quickly transferred to the disc stars. As a consequence a 
substantial fraction of disc is ejected into the stellar halo. In 
our simulations the stellar halo is defined as all the stellar 
content present outside three scale radii of the central disc 
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Figure 7. The fraction of the central stellar disc mass that ends 
up in the bulge of the galaxy as a function of the baryonic merger 
ratio. The blue points are our simulation results and the green 
line is the fit to our result. 

and/or 5 kpc above or below it. Figure [H] shows the fraction 
of central disc stellar mass that is dispersed into the stellar 
halo of the galaxy as a function of the baryonic merger ratio. 

A disc destroyed during a major merger does not en¬ 
tirely end up in the newly formed bulge but a significant 
fraction of its mass, up to ~ 22%, is ejected into the halo. 
If this mass transfer from the central disc to the halo is 
neglected then the bulge fraction of galaxie s will be signif¬ 
icantly overestimated. Recent results from iBernardi et al.l 
ll2013l) have shown that there is a substantial amount of 
starlight in the extended envelopes of massive galaxies. Our 
results show that mergers might be one of the ways to cre¬ 
ate massive amounts of intercluster light in these galaxies. 
Empirically we find that the fraction of central disc mass 
dispersed into the stellar halo during a merger (/d/i)is given 
by 

fdh = 0.22fib (13) 

These results show a revised picture of the the fate of 
a central stellar disc during a major merger. For a binary 
merger, about 40% of the disc mass loses angular momentum 
and is then transferred into the bulge. More than 20% the 
initial mass gains enough energy to escape from the central 
region and is dispersed into the galactic halo. Finally a small 
part of the disc is able to survive the merger and form a 
smaller (thicker) disc structure around the newly formed 
bulge. These results show that mergers are not as effective 
as previously thought in creating galactic bulges. 

3.5 Sites of star formation 

There has been a lot of effort put in to quantifying t he ef- 
ficiency of merger d riven starbursts (ICox et al.ll2008l . H09, 
iKarman et al.ll2015ll . H09 assume that all the disc gas mass 
that enters the bulge during a merger will be available for 
star formation and there is a p articular efficienc y for the con¬ 
version of this gas into stars. ICox et ^ ll2008ll have looked 
at star burst efficiency of the entire galaxy. The SAMs use 



Figure 8. The fraction of the central disc mass transferred to the 
halo as a function of the baryonic merger ratio. The green curve 
is out fit to the data. A significant fraction of the stars end up in 
the stellar halo of the central galaxy. 

these results for starburst efficiency but they assume that 
star formation mainly takes place in the bulge of the galaxy. 

Recent simulations have shown th at there is extende d 
clumpy star formation in major mergers (IPowell et al.ll20i^ i. 
They find that most mergers have an extended star forma¬ 
tion component during the early stages of a merger, but star 
forma tion becomes nuclear as the galaxies approach coales¬ 
cence. iMoreno et al.l ll2015l i find that the enhanced star for¬ 
mation in mergers is a combination of intense enhancements 
within the central kpc and moderately suppressed activity 
at larger galacto-centric radii. Here we try to distinguish and 
quantify the amount of star formation in the bulge (central) 
and disc (extended) components. To do this we look at the 
amount of new stars formed at the end of the simulation in 
the remnant compared to the new stars formed in the isola¬ 
tion run. This additional star formation is then due to just 
the merger. Our definition of the SF enhancement include 
stars formed in the satellite as well. A more sophisticated 
analysis of the merger driven st ar formation rates is pre¬ 
sented in a complementary work llKarman et al.l[2015ll . but 
the general trends are investigated here. 

The top panel of Fig shows the fractional change in 
star formation rate, i.e., the star formation in merger sim¬ 
ulation minus star formation in isolated galaxy divided by 
the star formation in isolated galaxy, as a function of the 
merger ratio for both bulge (red circles) and disc (blue stars) 
component. We see that in most cases the star formation is 
enhanced both in the bulge and the disc. For very minor 
mergers we see an increase in SF by about 10% in both the 
bulge and disc, whereas in major mergers the SF in the bulge 
can be enhanced as high as 80 times the original value. On 
the other hand there is relatively low enhancement (of about 
20 times) in the SF in the disc during a merger. 

We can then turn the previous plot around and look at, 
for a particular enhancement in SF, what fraction of it takes 
place in the bulge and disc. The bottom panel of figO shows 
the fraction of SF enhancement which takes place in the 
bulge and disc. In most mergers most of the star formation 
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Table 3. Table listing the best fit parameters of the functional form f{x) = ax 


Transfer Channel 


Ub 

l^c 

expi 1. 9r}) 

Ph 

exp(1.9ri) 

exp( 1.9r?) 

- 

a 


a 

x^ 

a 


a 


a 

X^ 

a 

X^ 

fab 

0.92 

1.36 

1.01 

1.06 

0.91 

1.58 

1.71 

1.95 

1.89 

1.69 

1.86 

1.77 

fdb 

0.31 

0.26 

0.37 

0.19 

0.30 

0.32 

0.54 

0.25 

0.58 

0.26 

0.56 

0.29 

fdh 

0.20 

0.10 

0.22 

0.11 

0.20 

0.12 

0.31 

0.14 

0.35 

0.13 

0.34 

0.14 


0.67 

0.77 

0.99 

0.29 

0.64 

0.86 

1.23 

0.86 

1.82 

0.43 

1.29 

0.85 


“ All the functional forms include a multiplicative factor of (1 — fqas) 



the hot halo plays in these cosmologically motivated merger 
simulations. We perform a new simulation for Tree 215240 
(chosen for its conveniently large merger ratio 0.77) without 
including the hot halo (both in the models of the primary 
and the satellite galaxies), but keeping all other components 
the same (stars, cold gas and dark matter). 

Figure [TU] shows the comparison between the two runs 
for various properties as a function of time. Without a hot 
halo (green curves) the primary galaxy becomes completely 
bulge dominated after the merger, as shown in the upper 
left panel. The B/T ratio increases from a value of about 
« 0.8 in the run with the hot halo (blue curves) to « 1.0 
when the hot gas reservoir is removed. This is due to the 
fact that the hot halo replenishes the disc of the galaxy, 
keeping it gas rich as shown on the bottom left panel (gas 
fraction of primary disc). This in turn helps the gas disc to 
survive the merger (eg. H09 and Eq. [T0]of this work). This 
can be seen in the top right panel of Fig. 1101 We can also 
see that the stellar mass transfer to the bulge from the disc 
remains unaffected and does not depend on the gas fraction 
of the disc (bottom right panel). The hot halo also provides a 
continuous supply of cold gas for SF which allows for highe r 
rate of disc regrowth after merger llDe Lucia et al.l I 2 OIII '). 
All together it is clear that the inclusion of a hot halo has 
several effects on the the final bulge to total ratio, since this 
hot gas reservoir is crucial for keeping the disc gas rich and 
for disc regrowth after a merger. Hence this component is 
important and should be taken into account in numeric al 
studies of mergers of galaxies (see also iMoster et aDl2012ll . 


Figure 9. Top: The fractional change in star formation in the cen¬ 
tral bulge (red circles) and the central disc (blue stars). Bottom: 
The fractional distribution of new stars formed due to merger. 


occurs in the bulge, but quite a few merger scenarios do show 
that the star formation can be triggered in discs as well. Our 
results seem to suggest that the relative contribution of disc 
starburst to the overall star formation enhancement is larger 
for minor mergers. In order to confirm and further quantify 
this effect higher resolution simulations are needed that go 
beyond the scope of this paper. 


3.6 Hot halo 


Previously, simulations which have included hot haloes, 
have looked into its effects in both minor mergers and 
majo r mergers in non-c osmological setting (iMoster et al.l 
I 2 OIII. iMoster et al.ll20l3 ) and on the star formation rates 
I Kim et al.l 2009l) ~It is then interesting to check what role 


3.7 Effect of resolution 

To test the effect of resolution we simulated Tree 215240 
(/i = 0.88) at a higher resolution, containing two times more 
particles than the fiducial run, i.e., 10® stellar particles in 
the final merger remnant. This also decreases the softening 
lengths by a factor of ~ 1.4. All the other simulation param¬ 
eters are the same as used in the fiducial run. Fig. [TT] shows 
the bulge to total ratio for both the fiducial and high res¬ 
olution simulations. Encouragingly, the different resolution 
runs match quite well. There are however, some subtle dif¬ 
ferences between the runs. The disc of the higher resolution 
run seems to be a bit more stable to perturbations than the 
fiducial run. This is evidenced more clearly after the merger 
takes place. The disc of the primary in the fiducial run is 
excited and there seems to be a larger perturbations on the 
low resolution disc as evidenced by the larger variation of 
the B/T ratio on short timescales. The general behaviour, 
however, is quite well converged. Although, we only show the 
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Figure 10. Effect of the hot halo on the BjT, fraction of satellite mass given to central bulge and fraction of central disc mass given to 
bulge. The blue curve indicates a simulation including a hot gaseous halo and the green curves are simulations without the hot halo. 


comparison of the evolution of the B/T ratio, we checked the 
other mass transfer channels and they all seem to be well 
converged with respect to resolution as well. 


4 DISCUSSION AND CONCLUSIONS 

Theoretical models of galaxy formation and evolution usu¬ 
ally assume that star formation in galaxies mainly occurs 
in a disc-like structure as a consequence of the conservation 
of angular momentum acquired by their host DM haloes. In 
such a scenario spheroidal structures like bulges are created 
by dissipative mechanisms, like galaxy mergers and disc in¬ 
stabilities, which remove angular momentum from stars and 
gas, funnelling them towards the centre of galaxy. Most the¬ 
oretical models, like SAMs, employ simplified formulae for 
dealing with these mass transfer processes, in the idealized 
approximation of a sequence of binary mergers. 

In this paper, we analyse in full detail the physics 
of galaxy mergers, and how mergers change the morphol¬ 
ogy of galaxies, with particular emphasis on the mass 
transfer processes between the different components of the 
merging galaxies. Our appr oach (first presented and de¬ 
scribed in iMoster et al.ll2014li consists of using high resolu¬ 
tion smoothed particle hydrodynamic simulations of galaxy 


mergers where the initial conditions are taken from a com¬ 
bination of cosmological realization of dark matter haloes 
merger trees and semi-analytic models. The cosmological 
merger trees ensure that the orbits and merger timings of 
DM substructures are cosmologically consistent, while the 
SAM modelling provides a reasonable guess for the proper¬ 
ties of galaxies living in these haloes. 

We run a series of simulations from z = 1 to z = 0, 
including both single and multiple mergers and we quantify 
the amount of morphological transformation in these model 
galaxies. We decompose the bulge and disc component of 
model galaxies using all the six dimensional phase space of 
position and velocity and considering the mass distribution 
as a function of the rotational support. This decomposition 
approach is very robust and allows us to track mass transfers 
from the satellite to the bulge, from the initial disc to the 
bulge, the mass of the disc dispersed into the stellar halo of 
the galaxy and the B/T evolution. This then allows us the 
quantify the amount of mass transfer in each channel as a 
function of the merger parameters. 

Mass transfers from the satellite to the bulge: In ma¬ 
jor mergers most of the satellite’s mass is deposited into the 
bulge of the central galaxy. On the other hand, a small satel¬ 
lite galaxy deposits most of its baryonic mass into the halo 
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Figure 11. The evolution of B/T ratio as a function of time for 
Tree 215240 {fi = 0.88) in both the low and high resolution run. 
The high resolution run has two times as many particles as the 
fiducial run and has ~ 1.4 times higher spatial resolution. 

(due to stellar stripping and ram pressure) before its final 
coalescence with the central object. 

Mass transfers from the central disc to the central bulge: 
The interaction with the incoming satellite causes matter 
inflows from the central disc to the central bulge. We confirm 
the findings of H09, that a gas rich disc is able to survive 
mergers more effectively. We note here that we do not resolve 
the small scale physics of t he ISM, instead choosing to model 
it using an effective EOS llSpringel &: Hernauistll^OSl ). This 
model provides a high thermal pressure support in the ISM, 
which smoothes out any instabilities in the disc. This model 
works quite well for low redshift redshift galaxies which are 
not highly star forming. At high redshifts the ISM is usually 
very turbulent and clumpy. Recent studies have shown that 
masses and sizes of disc components that survive a major 
merger are severely reduced when the hydr odynamics of the 
high redshift ISM are modelled properly llBournaud et al.l 
[m3). Since the starting point of our simulations is a = 1, 
the disc is expected to be much more calm. Therefore the 
ISM model used in our work is quite sufficient to capture 
the physics of mergers at the redshifts we consider. 

However, our results show that the amount of central 
stellar disc mass transferred to the central bulge is quite low. 
Even for 1 : 1 mergers only 37% of the central stellar disc is 
given to the central bulge. In addition a fraction of the cen¬ 
tral disc mass is expelled outwards due to close gravitational 
encounters (about 22% for 1 : 1 mergers), adding up to the 
diffuse stellar halo population. This mass transfer channel is 
generally poorly modelled in many SAMs. 

Merger driven star formation: Our results also confirm 
that there is a sharp enhancement of star formation during 
a merger. Many SAMs assume that all of the merger driven 
starburst takes plac es in the bulge. However, our simulation s 
echo recent results (|Powell et ai.ll201^ : iMoreno et al.ll201!Th . 
which show that there is an enhancement of star formation in 
the discs of central galaxies in addition to the enhancement 
in the bulge. As expected, larger starbursts are attributed 
to major mergers, but the relative contribution of disc star- 
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burst to the overall star formation enhancement is larger for 
minor mergers. 

Effect of hot halo: One of the major improvements of 
this work, with respect to previous similar attempts, lies in 
the modelling of a hot gaseous halo component: the existence 
of such a component (in addition to the cold phase present 
in the disc of the galaxy) is predict e d by cosmological sim- 


ulations fe.e. Sommer-LarserJ 

20061: Johansson et alJ l2009l: 

IPasmussen & Ponmanll2009l: ll 

ansen et al.ll201l|) and obser- 

vations le.g. Owen & Warwick 

20091. Anderson et al.ll2013l). 


The hot gaseous halo affects the final bulge to total ratio 
of the galaxy: in fact it is able to replenish the galaxy disc 
with fresh infalling cold gas, keeping it gas rich and fuelling 
disc star formation over longer time scales. 

We thus propose a series of new fitting formulae able to 
capture the trends in the relations between the amplitude 
of mass transfers and baryonic (stars -|- cold gas) merger 
ratio. We are well aware of the limitation of our analysis, a 
relatively small sample of simulated haloes. More work is of 
course needed to increase the simulation sample, including 
more merger trees at different mass scales and widening the 
range of orbital/merger parameter. This larger sample would 
then provide us with better constrains on both the mean 
relations and their scatter, as well as a testbed for extreme 
cases, such as multiple merger scenarios. Despite the small 
sample considered in this paper, our results are indicative 
of a need for revision in our understanding of mass flows 
involved in a galaxy merger. In a work (Fontanot et ah, in 
preparation), we plan to include our fitting formulae in state- 
of-the-art SAMs, to study the implications of our findings on 
a cosmologically significant galaxy sample, in determining 
the fraction of bulge and disc dominated galaxies. 

Overall our study shows that the combined effect of 
satellite accretion and disc destruction in building bulges is 
not as efficient as previously thought and this will possibly 
alleviate some of the tensions between the observed frac¬ 
tion of bulgeless galaxies and the hierarchical scenario for 
structure formation as predicted by the ACDM model. 
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